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Abstract 

This paper reports the tenth-order QED contribution to the lepton g—2 from the gauge-invariant 
set, called Set III(c), which consists of 390 Feynman vertex diagrams containing an internal fourth- 
order light-by-light-scattering subdiagram. The mass-independent contribution of Set ni(c) to the 
electron g — 2 (og) is 4.9210 (103) in units of (a/vr)^. The mass-dependent contributions to 
from diagrams containing a muon loop is 0.00370 (37) (q/tt)^. The tan- lepton loop contribution 
is negligible at present. Altogether the contribution of Set ni(c) to is 4.9247 (104) (a/vr)^. We 
have also evaluated the contribution of the closed electron loop to the muon g—2 (a^). The result 
is 7.435 (134) {a/irf. The contribution of the tau-lepton loop to is 0.1999 (28) {a/irf. The 
total contribution of various leptonic loops (electron, muon, and tau-lepton) of Set III(c) to is 
12.556 (135) {a/irf. 

PACS numbers: 13.40.Em, 06.20.Jr,12.20.Ds,14.60.Cd 
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I. INTRODUCTION 



The anomalous magnetic moment = {g—2)/2 of the electron has played the central role 
in testing the validity of quantum electrodynamics (QED) as well as the Standard Model. 



On the experimental side, the 



the precision of 0.24 x 10 



atest measurement of Og by the Harvard group has reached 
3: 



ae(HV08) = 1 159 652 180.73 (0.28) x 10"^^ [0.24ppb] . 



(1) 



T 

B 



le theoretical prediction thus far consists of QED corrections of up to the eighth order 



-|5|, direct evaluation of hadronic corrections [6 
down from their contributions to the muon g — 2 



-Il2|. and electroweak corrections scaled 
13l-|l5|. To compare the theory with the 



measurement ([1]), we also need the value of the fine structure constant a determined by 
a method independent of (7 — 2 . The best value of such an a available at present is one 
obtained from the measurement of h/rriYib, the ratio of the Planck constant and the mass 
of Rb atom, combined with the very precisely known values of the Rydberg constant and 
WRb/me: [l6| 



a-^(RblO) = 137.035 999 037 (91) [0.66ppb]. 



With this a the theoretical prediction of a^. becomes 



ae(theory) = 1 159 652 181.13 (0.11)(0.37)(0.02)(0.77) x 10 



(2) 



-12 



(3) 



where the first, second, third, and fourth uncertainties come from the calculated eighth- 
o.d. QED te.. a, the cude .e,.h-o.de. e.i.a.e Q. .he had™.c and e.eCwea. 
contributions, and the fine structure constant ([2]), respectively. The theory (|3]) is in good 
agreement with the experiment (II]): 



ae(HV08) - ae(theory) = -0.40 (0.88) x 10 



-12 



(4) 



proving that QED (Standard Model) is in good shape even at this very high precision. 

Eq. OH]) shows clearly that the largest source of uncertainty is the fine structure constant 
02]). To put it differently, a non-QED a, even the best one available at present, is too 
crude to test QED to the extent achieved by the theory and measurement of Og. Thus it 



2 



makes more sense to test QED by an alternate approach, namely, obtain a from theory and 
measurement of a^^: 



a 



-i(a,08) = 137.035 999 085 (12)(37)(2)(33) [0.37ppb], (5) 



where the first, second, third, and fourth uncertainties come from the calculated eighth-order 
QED term, the crude tenth-order estimate, the hadronic and electroweak contributions, and 
the measurement of ae(HV08), respectively. 

Although the uncertainty of a~^(ae08) in (|5]) is a factor 2 smaller than a~^(RblO), it is 
not a firm factor since it depends on the estimate of the tenth-order term, which is only a 
crude guess |l7|. For a more stringent test of QED, it is obviously necessary to calculate 
the actual value of the tenth-order term. To meet this challenge we launched several years 



18|-|20| 



ago a systematic program to evaluate the complete tenth-order term 

The lOth-order QED contribution to the anomalous magnetic moment of an electron can 
be written as 



-TT 



where mjm^ = 4.836 331 66 (12) x 10"^ and m^/m^ = 2.875 64 (47) x 10"^ tl7|]. In the 
rest of this article the factor (^j is suppressed for simplicity. 

The diagrams contributing to the mass-independent term A^^^^ can be classified into 
six gauge-invariant sets, further divided into 32 gauge-invariant subsets depending on the 
nature of closed lepton loop subdiagrams. Thus far, numerical results of 30 gauge-invariant 
subsets, which consist of 5928 vertex diagrams, have been published 



21 l2lH27[^r submitted 



for publication [28| . Five of these 30 subsets are also known analytically |29|, l30| . They are 
in good agreement with our calculations. 

In this paper we report the contribution to from the gauge-invariant subset called 
Set III(c), which consists of 390 vertex diagrams and is represented by 24 self-energy- like 
diagrams of Figure [TJ A characteristic feature of these diagrams is that they have vertex 
and self-energy subdiagrams which contain a light-by-light (/-/) scattering subdiagram. They 
can be classified into three types: 

1. Sixth-order vertex subdiagrams containing an /-/ loop externally as is shown in Fig- 
ure [21 Here "external" means that one of the photons is external to the subdiagrams 
containing the /-/ loop. 




A5 
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FIG. 1: Tenth-order self-energy- like diagrams in which lepton lines propagate in the magnetic field. 
They represent 390 vertex diagrams of Set III(c). Assignment of Feynman parameters zi, Z2, ■ ■ ■ , 
to the lepton lines and Za-, Zb, .., Ze to the photon lines is indicated in the figure Al, Bl, and CI. 

2. Eighth-order vertex subdiagrams containing an /-/ loop externally as is shown in Fig- 
ure [31 They are obtained by applying a virtual photon correction on the open fermion 
line of a sixth-order /-/ vertex subdiagram of Figure [2J 



3. Eighth-order vertex and self-energy subdiagrams which contain an /-/ loop internally. 
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FIG. 2: Vertex diagrams of sixth order containing an l-l subdiagram. The diagram 6LL(3) and 
6LL(4) are identical with each other because of the time-reversal symmetry. There are six vertex 
diagrams of this type, taking into account of two directions in which the closed fermion loop can 
take. 

This type appears for the first time in the tenth-order perturbation theory of QED. 
See Figure ID 

The vertex subdiagrams of type 1 and type 2 do not have their Ward-Takahashi-related 
self-energy subdiagrams which vanish identically due to Furry's theorem. Thus, the gauge- 
invariant sums of the vertex renormalization constants of these external /-/ subdiagrams also 
vanish identically. 

For vertex subdiagrams of type 3, the corresponding self-energy subdiagrams do exist 
which have an internal /-/ diagram. In this case both vertex subdiagram and self-energy 
subdiagram have UV divergence due to the sixth-order external /-/ vertex subdiagram. 
Because of these specific features of an /-/ scattering diagram and vertex diagrams Con- 
or the Set III(c) an approach different from the one 
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20| . Our formulation and treatment of UV diver- 



taining an external /-/ loop, we adopt 
used for a diagram without an /-/ loop 
gences and IR divergences due to subdiagrams are described in Sec. [Tll Results of numerical 
evaluation will be presented in Sec. IIIII Sec. HV] is devoted to the summary and discussion 
of this work. Renormalization of these diagrams is described in Appendix |A1 

II. FORMULATION 

Instead of dealing with the 390 vertex diagrams of Set III(c) individually, we consider 
the sum A'^ of a set of vertex diagrams that are obtained from a self-energy-like diagram 
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FIG. 3: Various diagrams of the eighth order needed for renormalization. 8LLa (a = E, F, G, H, I) 



is denoted as LLa in Ref. 



. A vertex diagram 8LLa{i) {i = 5, 6, 7) is obtained by inser ting an 



external vertex into a fermion Hne i of the diagram 8LLa. It is denoted as LLa{i) in Ref. 



32]. 
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FIG. 4: The eighth-order self-energy-hke diagrams 8LLJ, 8LLK, and 8LLL, containing a hght- 
by-light scattering loop internally. Open lepton lines propagate in the weak magnetic field. They 
represent 18 vertex diagrams in total. The vertex diagram 8LLJ{1) is a part of the diagram 8LLJ 
in which the magnetic vertex is attached only to the fermion line 1 of the self-energy diagram 
8LLJ. 

of Figured] by inserting a magnetic vertex j'^ in the lepton lines 1, 2, 3, 4, and 5. We 
rewrite this as 



dA^{p,q) 
dqu 



q=0 



dpu 



(7) 



with the help of the Ward-Takahashi identity, where p—q/2 and p+q/2 are the 4-momenta of 
incoming and outgoing lepton lines and (p — = + = m^. Each sum corresponds 
to one of the 24 self-energy-like diagrams shown in Figured] The g—2 term is projected out 
from the right-hand side of (^^. 

A. Construction of Unrenormalized Integrals 



Each diagram Q of Figure d] can be expressed by a momentum integral applying the 
Feynman-Dyson rule. Introducing Feynman parameters zi, Z2, ■ ■ ■ , Zg for the electron prop- 
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agators and Za, Zb, . . . , Ze for the photon propagators (see the figures Ai, Bi, and Ci of 
Figured]), we carry out the momentum integration analytically by means of a home-made 
program written in FORM 33|. This leads to an integral of the form 

Mg = -(^)\\ [{dz)g 



4 



En 



Co ^ El + Ci ^ 



No + Zo Ni + Zi 



+ 



where En, C„, Nn and Z„ are functions of Feynman parameters, and "symbolic" building 



blocks Ai, Bij, Cij, for i, j = 1,2, . . . ,9. n is the number of contractions (see 3J] for defini- 



tions). U is the Jacobian of transformation from the momentum space variables to Feynman 
parameters. Ai is the scalar current defined by 



Ai 



Vi 



u 



1 for i 



1,2,3,4,5 



for i = 6, 7, 8, 9 



(9) 



and 



{dz)g = l[dz,{l-Y.^^y (10) 

See, for example, Ref. [19| for definitions of Bij and Cij. V is obtained by combining all 
denominators of propagators into one with the help of Feynman parameters. It has a form 
common to all diagrams of Figure [1] 

9 e 

V = Y,z,{l-Ai)ml + Y.'^^l^ (11) 

j=l k=a 

where rui and are the rest masses of lepton i and photon k, respectively. Of course, rui 
is independent of i and is independent of k. But it is useful to distinguish different 
lepton lines and photon lines in deriving Eq. ([H]). The form of Ai as a function of Feynman 
parameters depends on the structure of individual diagram Q of Figure [H 



B. Renormalization 



The diagrams of Set III(c) as a whole form a (formal) gauge-invariant set. However, 
individual diagrams have UV divergences arising from the light-by-light-scattering (/-/) sub- 
diagram as well as vertex subdiagrams or self-energy subdiagrams. All these divergences 
must be regularized in advance. In order to maintain gauge invariance the /-/ subdiagram 
may be regularized by the Pauli-Villars method or by the dimensional regularization. The 
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sixth-order vertex renormalization constant associated with the diagram containing an 
subdiagram and the eighth-order vertex renormahzation constant containing an /-/ sub- 
diagram are logarithmically divergent, but their sum over all diagrams vanishes due to 
Ward-Takahashi identity. Note that the self-energy diagrams associated with these vertex 
diagrams do not exist in QED because of the Furry's theorem. 

As is indicated in the figures Al, Bl, CI of Figured! we denote open fermion lines as 1, 2, 
3, 4, 5, fermion lines forming a closed loop as 6, 7, 8, 9, and photon lines as a, b, c, d, e. We 
will identify a subdiagram containing open lepton lines in terms of a subset of (1,2,3,4,5). 
For instance, the vertex subdiagram (1,2) of Al will be denoted by (1,2), and the vertex 
subdiagram {2,3,4,5; 6,7,8,9; b,c,d} of Al will be denoted by (2,3,4,5). An exception is the /-/ 
subdiagram, which will be denoted as (6,7,8,9). Under this convention the diagram Al has 
five divergent subdiagrams (1, 2), (4, 5), (1, 2, 3, 4), (2, 3, 4, 5), and (6, 7, 8, 9). The fifteen UV 
subtraction terms can be constructed from these subdiagrams following the Zimmermann's 



forest formula 35 1. 



Diagrammatically, the second-order vertex subdiagram appears not only in the forests 
including the subdiagram (1,2) but also in the forest (2, 3, 4, 5)(4, 5). In the latter, the 
reduced diagram (2, 3) forms a second-order vertex diagram. We will treat renormalization 
of this implicit second-order vertex in a manner different from the explicit second-order 
vertex. A detailed account will be given in Appendix A. 

The UV divergence arising from the explicit second-order vertex (1,2) of the diagram Al 
can be subtracted by an integral defined by the i<'i2-operation [3J] applied on the integral 
Mm- The iri2-operation is defined in such a way that the result of the operation factorizes 
exactly as 

i^,2M^, = 4^'^)^^M£tj\ (12) 

where L^''^^''^^ is the UV-divergent part of the second-order on-shell vertex renormalization 
constant L^''^^ and Mg^'^j^ is the magnetic moment amplitude from the eighth-order self- 
energy-like diagram 8LLJ of Figure HI 

UV divergences from the explicit second-order vertex subdiagram are also found in the 
diagrams B1,C1, A6, B6, and C6. UV divergences due to the explicit second-order self- 
energy-like subdiagram come from the diagrams ^45, B5, C5, A8, B8, and C8. The renor- 
malization scheme in which only these second-order divergences appear are handled by the 
i^-operation and is described in Appendix \M 
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All other subdiagrams contain an subdiagram, which we treat by the Pauli-Villars 
method or by the dimensional regularization. For instance, in the latter method, let Fai{d) be 
one of such integrals defined in d dimension, where ai takes values a = A,B,C; i = 1,2,..., 8. 
Let Gai{d) be Fai{d) in which the /-/subdiagram (of the form Il^,^„p{ki, ^2, ^3, ^4)) is replaced 
by the tensor with zero external momenta, namely, n^;^o-p(0, 0, 0, 0). Let us rewrite Fai{d) 
symbolically as 

[F^^{d)-GUd)] + GUd), (13) 

where by "symbolically" we mean that subtraction is performed on the integrand before the 
integration is carried out. Now we can safely take the limit d ^ A for the term [Fai{d) — 
Gai{d)] since its integrand does not cause UV divergence. Of course, the second term Gai{d) 
is singular for d — t- 4. However, gauge invariance guarantees that the sum of Gai{d) over all 
diagrams of Figure [1] vanishes for any value of dimension d: 

C 8 

$^$^r/,G„.(ci) = 0, (14) 

a=A i=l 

where rji = 2 ior i = 4, 7, 8, and r/j = 4 for i = 1, 2, 3, 5, 6. Thus, in the end, we have to 
compute only 

\im[F^,{d)-G^,{d)]. (15) 

Of course the same result is obtained by the Pauli-Villars method. To avoid crowded nota- 
tions let us use Fai{4:) instead of Eq. f|T5|) in the following. 

Each self-energy-like diagram of Figure [1] represents the sum of five vertex diagrams. 
Diagrams obtained by reversing the momentum fiow within the /-/ loop are not shown but 
they give the same integrals as the original ones. Another factor 2 must be included for 
diagrams that are not symmetric under time-reversal. Thus, integrals for diagrams such as 
Al actually represent 2x2x5 vertex diagrams. The g — 2 contribution from the sum of 
all diagrams of Set III(c), after the renormalization described in Appendix |X] is carried out, 
can thus be written as 

c 8 

4'°) [Set III(c)('^'^)] = ^ ^ry,AM(f ^) - SALSaAM^M, (16) 

a=A i=l 

where /i refers to the open lepton line and I2 refers to the closed lepton line, rji = 2 for 
z = 4, 7, 8, and r]i = 4 for i = 1, 2, 3, 5, 6. ALB2 and AM^ji^l are defined in Appendix [Al 
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III. NUMERICAL RESULTS 



Evaluation of integral AM^'/'^^ is carried out by the adaptive- iterative Monte-Carlo in- 



tegration routine VEGAS 



The results for the case = (ee) are listed in Table [H 



From this table and Table HIl listing the residual renormalization terms we obtain 

AS'°^[Set III(c)(''^^] = 4.9210 (103). (17) 

The contribution of the muon loop to Oe can be calculated from the data listed in Table 
inland Table [TTl 

4'°^ [Set III(c)('^"^)] = 0.00370 (37). (18) 

The contribution of the tau-lepton loop to Og is within the uncertainty of fll7p . Thus the 
total QED contribution to ai^°^ is essentially the sum of flTTl) and flTSl) : 
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i(^°)[Set 111(c)] = 4.9247 (104) . (19) 



FORTRAN programs for can be readily adapted to the evaluation of a^. The results 
of evaluation of the contribution of the electron loop to the muon g — 2 are listed in Table 
IIV[ From this table and Table [TTl we obtain 

4'°^ [Set lll(c)('"')] = 7.435 (134). (20) 

The contribution of the tau-lepton loop to is calculated from the data listed in Table El 
and Table ITTl 

4'°^ [Set 111(c) = 0.1999 (28). (21) 
The total QED contribution to 4^°^ is the sum of ffTT]) . f l20|) . and (|2T]): 

a(,i°)[Set 111(c)] = 12.556 (135) . (22) 



IV. DISCUSSION 

All programs of diagrams of the Set 111(c) were written in two independent ways, in order 
to detect possible programming error. No such error was found. 

The value of A2^'''*[Set III(c)'^™'^'*] given in (120!) is not much larger than that of 
^S'°^[Set III(c)('^')] given in ( ITTj) . This is somewhat unexpected since, as is seen from Table 
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TABLE I: Contributions of diagrams of Set in(c) to Og for (^1/2) = (ee). The superscript (ee) is 
suppressed for simplicity. The multiphcity ni? is the number of vertex diagrams represented by the 
integral and is incorporated in the numerical value. All integrals are evaluated initially with 10^ 
sampling points per iteration, iterated 50 times, followed by 10^ points, iterated several times. 



Integral 


np 


Value (Error) 
including np 


Sampling per 
iteration 


No. of 
iterations 


AMai 


20 


-4.255 92 


(253) 


1 X 10^, 


1 X 10^ 


50, 400 


AMa2 


20 


4.938 78 


(244) 


1 X 10^, 


1 X 10^ 


50, 300 


AMa3 


20 


-1.546 88 


(246) 


1 X 10^, 


1 X 10^ 


50, 345 


AMa4 


10 


-0.323 88 


(127) 


1 X 10^, 


1 X 10^ 


50, 30 


AMa5 


20 


6.320 29 


(153) 


1 X 10^, 


1 X 10^ 


50, 60 


AMa6 


20 


-5.660 33 


(218) 


1 X 10^, 


1 X 10^ 


50, 300 


AMa7 


10 


2.284 61 


(173) 


1 X 10^, 


1 X 10^ 


50, 65 


AMa8 


10 


1.362 06 


(129) 


1 X 10^, 


1 X 10*^ 


50, 20 


AMbi 


20 


5.693 53 


(293) 


1 X 10*^, 


1 X 10^ 


50, 412 


AMb2 


20 


-7.018 17 


(273) 


1 X 10^, 


1 X 10^ 


50, 302 


AMb3 


20 


3.735 46 


(260) 


1 X 10^, 


1 X 10^ 


50, 342 


AMba 


10 


-0.052 76 


(122) 


1 X 10^, 


1 X 10^ 


50, 30 


AMb5 


20 


-4.739 40 


(166) 


1 X 10^, 


1 X 10*^ 


50, 60 


AMbg 


20 


3.061 01 


(212) 


1 X 10^, 


1 X 10^ 


50, 300 


AMb7 


10 


0.351 39 


(168) 


1 X 10^, 


1 X 10^ 


50, 65 


AMb8 


10 


-0.793 52 


(136) 


1 X 10^, 


1 X 10^ 


50, 20 


AMci 


20 


0.377 40 


(279) 


1 X 10^, 


1 X 10*^ 


50, 417 


AMc2 


20 


3.054 41 


(241) 


1 X 10^, 


1 X 10*^ 


50, 300 


AMc3 


20 


-1.329 04 


(260) 


1 X 10^, 


1 X 10^ 


50, 338 


AMcA 


10 


0.435 88 


(131) 


1 X 10^, 


1 X 10^ 


50, 30 


AMc5 


20 


-3.729 22 


(159) 


1 X 10^, 


1 X 10^ 


50, 60 


AMc6 


20 


4.273 41 


(258) 


1 X 10^, 


1 X lO'^ 


50, 300 


AMci 


10 


-2.233 00 


(159) 


1 X 10^, 


1 X 10*^ 


50, 65 


AMc8 


10 


-1.514 28 


(142) 


1 X 10^, 


1 X 10*^ 


50, 20 



TABLE II: Auxiliary integrals for Set III(c). Some integrals are known exactly. Other integrals 
are obtained by the integration routine VEGAS. The superscript (/i I2) indicates that the open 
and closed fermion lines consist of fermions l\ and Z2, respectively. The letters e, m, and t stand 
for electron, muon, and tau-lepton, respectively. 



Integral 



Value (error) 



Integral 



Value (error) 



M2 



0.5 



ALBo 



0.75 



AM 



(ee) 



8JKL 



AM 



(em) 



8JKL 



-0.990 72 (11) 
-0.000 177 8 (13) 



AM 



(me) 



%JKL 



AM 



(mt) 



8JKL 



-4.432 43 (59) 
-0.015 87 (5) 
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TABLE III: Contributions of diagrams of Set III(c) to ae for (/1/2) = (em). The superscript (em) is 
suppressed for simplicity. The multipHcity n^? is the number of vertex diagrams represented by the 
integral and is incorporated in the numerical value. All integrals are evaluated with 10^ sampling 
points per iteration, iterated 50 times, and subsequently evaluated with 10^ sampling points per 
iteration, iterated 50 times. 



Integral 


np 


Value (Error) 
including np 


Sampling per 
iteration 


No. of 
iterations 


AMai 


20 


-0.016 78 


(15) 


1 X lO'^, 


1 X 10* 


50, 50 


AMa2 


20 


-0.004 71 


(8) 


1 X 10^, 


1 X 10*^ 


50, 50 


AMa3 


20 


0.000 99 


(6) 


1 X 10^, 


1 X 10* 


50, 50 


AMa4 


10 


-0.003 93 


(1) 


1 X 10^, 


1 X 10* 


50, 50 


AMa5 


20 


0.007 01 


(1) 


1 X lO'^, 


1 X 10* 


50, 50 


AMa6 


20 


-0.023 43 


(12) 


1 X 10^, 


1 X 10* 


50, 50 


AMa7 


10 


-0.001 00 


(2) 


1 X 10^, 


1 X 10* 


50, 50 


AMa8 


10 


0.001 97 


(1) 


1 X 10^, 


1 X 10* 


50, 50 


AMbi 


20 


0.007 61 


(15) 


1 X 10^, 


1 X 10* 


50, 50 


AMb2 


20 


0.000 37 


(8) 


1 X lO'^, 


1 X 10* 


50, 50 


AMb3 


20 


0.000 46 


(4) 


1 X lO'^, 


1 X 10* 


50, 50 


AMb4 


10 


0.003 05 


(1) 


1 X 10^, 


1 X 10* 


50, 50 


AMb5 


20 


0.010 68 


(1) 


1 X 10^, 


1 X 10* 


50, 50 


AMb6 


20 


0.015 17 


(11) 


1 X 10^, 


1 X 10* 


50, 50 


AMb7 


10 


0.002 24 


(2) 


1 X lO'^, 


1 X 10* 


50, 50 


AMb8 


10 


-0.013 72 


(1) 


1 X lO'^, 


1 X 10* 


50, 50 


AMci 


20 


0.010 57 


(12) 


1 X 10^, 


1 X 10* 


50, 50 


AMc2 


20 


0.004 88 


(5) 


1 X 10^, 


1 X 10* 


50, 50 


AMc3 


20 


0.000 87 


(4) 


1 X 10^, 


1 X 10* 


50, 50 


AMc4 


10 


0.000 84 


(1) 


1 X lO'^, 


1 X 10* 


50, 50 


AMc5 


20 


-0.018 16 


(1) 


1 X 10^, 


1 X 10* 


50, 50 


AMc6 


20 


0.009 57 


(10) 


1 X 10^, 


1 X 10* 


50, 50 


AMc7 


10 


0.000 93 


(2) 


1 X 10^, 


1 X 10* 


50, 50 


AMcs 


10 


0.011 40 


(1) 


1 X 10^, 


1 X 10* 


50, 50 



IIVI individual integrals contributing to A2 [Set 111(c) ] are an order of magnitude larger 
than those given in Table HI Presumably, the modest value of ( 120|) is a consequence of strong 
cancellation among contributing integrals. 
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TABLE IV: Contributions of diagrams of Set III(c) to for (/1/2) = (me). The superscript (me) 
is suppressed for simplicity. The multiphcity is the number of vertex diagrams represented by 
the integral and is incorporated in the numerical value. All integrals are evaluated initially with 
10^ sampling points per iteration, iterated 50 times, followed by 10^ points, iterated several times. 



Integral 


np 


Value (Error) 
including np 


Sampling per 
iteration 


No. of 
ixeraxions 


AMai 


20 


-18.722 


(37) 


1 X 10^, 


1 X 10*^ 


50, 376 


AMa2 


20 


40.155 


(26) 


1 X 10^, 


1 X 10^ 


50, 300 


AMa3 


20 


-3.780 


(36) 


1 X 10^, 


1 X 10^ 


50, 366 


AMaa 


10 


-18.309 


(11) 


1 X 10^, 


1 X 10^ 


50, 30 


AMa5 


20 


9.416 


(14) 


1 X 10^, 


1 X 10^ 


50, 80 


AMaq 


20 


-37.911 


(30) 


1 X 10^, 


1 X 10*^ 


50, 301 


AMa7 


10 


19.431 


(16) 


1 X 10®, 


1 X 10^ 


50, 85 


AMas 


10 


10.371 


(7) 


1 X 10^, 


1 X 10^ 


50, 70 


AMbi 


20 


54.402 


(38) 


1 X 10®, 


1 X 10^ 


50, 471 


AMb2 


20 


-73.374 


(29) 


1 X 10®, 


1 X 10^ 


50, 300 


AMb3 


20 


29.954 


(38) 


1 X 10®, 


1 X 10*^ 


50, 382 


AMba 


10 


2.578 


(13) 


1 X 10®, 


1 X 10^ 


50, 30 


AMb5 


20 


-49.408 


(20) 


1 X 10®, 


1 X 10^ 


50, 80 


AMb6 


20 


-1.509 


(33) 


1 X 10®, 


1 X 10^ 


50, 301 


AMbi 


10 


9.521 


(20) 


1 X 10®, 


1 X 10^ 


50, 85 


AMb8 


10 


29.116 


(9) 


1 X 10®, 


1 X 10*^ 


50, 70 


AMci 


20 


-31.212 


(37) 


1 X 10®, 


1 X 10^ 


50, 497 


AMc2 


20 


36.233 


(32) 


1 X 10®, 


1 X 10^ 


50, 300 


AMc3 


20 


-25.285 


(37) 


1 X 10®, 


1 X 10^ 


50, 409 


AMcA 


10 


15.428 


(16) 


1 X 10®, 


1 X 10^ 


50, 30 


AMc5 


20 


28.857 


(22) 


1 X 10®, 


1 X 10*^ 


50, 80 


AMcQ 


20 


43.793 


(38) 


1 X 10®, 


1 X 10*^ 


50, 310 


AMc7 


10 


-27.637 


(17) 


1 X 10®, 


1 X 10^ 


50, 85 


AMc% 


10 


-44.647 


(11) 


1 X 10®, 


1 X 10^ 


50, 75 
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RIKEN Integrated Cluster of Clusters (RICC) supercomputing systems, and the tp computer 
of Kobayashi-Maskawa Institute. 
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TABLE V: Contributions of diagrams of Set ni(c) to for (/1/2) = (mt). The superscript (mt) 
is suppressed for simplicity. The multiphcity np is the number of vertex diagrams represented by 
the integral and is incorporated in the numerical value. All integrals are evaluated initially with 
10^ sampling points per iteration, iterated 50 times, followed by 10^ points, iterated several times. 



Integral 


np 


Value (Error) 


Sampling per 


_ No. of 






including np 


iteration 


itcrcitions 


AMai 


20 


-0.423 43 (111) 


1 X 10®, 


1 X 10^ 


50, 40 


AMa2 


20 


-0.001 66 (64) 


1 X 10®, 


1 X 10^ 


50, 30 


AMa3 


20 


0.033 30 (59) 


1 X 10®, 


1 X 10^ 


50, 30 


AMa4 


10 


-0.102 91 (13) 


1 X 


10® 


50 


AMa5 


20 


0.327 57 (14) 


1 X 10®, 


1 X 10^ 


50, 15 


AMa6 


20 


-0.600 58 (75) 


1 X 10®, 


1 X 10^ 


50, 40 


AMa7 


10 


0.011 26 (30) 


1 X 10®, 


1 X 10^ 


50, 15 


AMas 


10 


0.065 01 (9) 


1 X 10®, 


1 X 10^ 


50, 15 


AMbi 


20 


0.246 29 (111) 


1 X 10®, 


1 X 10^ 


50, 40 


AMb2 


20 


-0.104 08 (64) 


1 X 10®, 


1 X 10^ 


50, 30 


AMb3 


20 


0.052 89 (49) 


1 X 10®, 


1 X 10^ 


50, 30 


AMba 


10 


0.076 46 (11) 


1 V 
1 A 


lU 


50 


AMb5 


20 


0.114 81 (13) 


1 X 10®, 


1 X 10^ 


50, 15 


A fl/f 

iaMbg 


on 


U.oy ( o\) \p( ) 


1 X 10®, 


1 X 10^ 


OU, 4U 


AMbi 


10 


0.043 30 (25) 


1 X 10®, 


1 X 10^ 


50, 15 


AMb8 


10 


-0.298 83 (9) 


1 X 10®, 


1 X 10^ 


50, 15 


AMci 


20 


0.251 61 (108) 


1 X 10®, 


1 X 10^ 


50, 40 


AMc2 


20 


0.137 58 (51) 


1 X 10®, 


1 X 10^ 


50, 30 


AMc3 


20 


-0.052 89 (51) 


1 X 10®, 


1 X 10^ 


50, 30 


AMci 


10 


0.026 64 (10) 


1 X 


10® 


50 


AMc5 


20 


-0.476 30 (12) 


1 X 10®, 


1 X 10^ 


50, 15 


AMc6 


20 


0.268 81 (65) 


1 X 10®, 


1 X 10^ 


50, 40 


AMc7 


10 


-0.038 53 (25) 


1 X 10®, 


1 X 10^ 


50, 15 


AMcs 


10 


0.213 27 (14) 


1 X 


10® 


50 


Appendix A: Renormalization of diagrams of Set III(c) 






Diagrams 


of Set III(c) 


, shown in Figure [H contain an 


/-/ subdiai 


^ram internally. Thus we 



find it convenient to pursue a renormalization scheme somewhat different from all other Sets 
contributing to the tenth-order g—2. As is indicated by figures Al, 51, CI of Figure [T], we 



denote open fermion lines as 1, 2, 3, 4, 5, closed fermion lines as 6, 7, 8, 9, and photon lines 
as a, b, c, d, e. We will identify a subdiagram containing open lepton lines in terms of their 
line numbers. For instance, the second-order vertex subdiagram {l,2;e} and sixth-order 
vertex subdiagram {4,5;6,7,8,9;b,c,d} of Al will be denoted by the superscript (1,2) and 
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(4,5), respectively. An exception is the subdiagram, which will be denoted as (6,7,8,9). 
Of course this is just for the sake of keeping track of where a particular subdiagram is 
located. The superscript will be removed when it is no longer needed. 



1. Al, Bl, CI 

Let us begin with the g—2 amplitude Mai- Noting that, out of 15 forests of the diagram Al 
mentioned in Sec. IIIBl 8 are hidden in our convention leading to Eq. (ITB]1 . the renormalized 
amplitude aAi can be written as 

_ M r (1,2) ,^(3,4,5) .(4,5) »^(1,2,3) . (1,2,3,4) ,.(5) . (2,3,4,5) ,.(1) 

+ L^"'4S,.,Mf ' + Ui^LSl^,Mi> + L?^'L^„,M^K (Al) 

As was discussed in Sec. [Tll all terms of ( lAlj) containing an subdiagram are to be under- 
stood as shorthands for the regularized quantity defined by Eq. ( ITSl) . In other words, the 
UV divergence arising from the /-/ subdiagram has been removed by the procedure described 
in Sec. |Tl] so that it can be treated as a UV-finite quantity. Msllj is the proper magnetic 



moment amplitude of the eighth-order diagram 8LLJ of Figure HI See [3l|, l32| for its precise 
definition. L2 is the vertex renormalization constant of the second order. L6LL(5) is the 
renormalization constant associated with the sixth-order vertex diagram 6LL{5) shown in 
Figure [2J Lsllf{7) and Lsllj{i) are the eighth-order vertex renormalization constants associ- 
ated with the self-energy-like diagrams 8LLF of Figure Eland 8LLJ of FigureHJ respectively. 

In the amplitude Mai, the ii'-operation is applied only on the explicit second-order vertex 
subdiagram (1,2). For other terms the full bodies of the vertex renormalization constants of 
the sixth- and eighth-orders are used and subtracted. These vertex renormalization constants 
are extracted from a vertex diagram Ti^{p, q), where {p — g/2)^ = (p + Q'/2)^ = rn'^, using the 
projection operator 

L = lTr[(^ + m)p'^r,]|,=o- (A2) 



The result is combined with the lower-order magnetic moment amp 



itude using, for in- 
stance, the factorization procedure described in Sec. Ill D of Ref. [19| backwards so that 
the combined formula is described by the same set of Feynman parameters as those of the 
unrenormalized magnetic moment Mai- Then the UV-finite quantity AM41 can be written 
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as 



- A/f rC/V(l,2) ,.(3,4,5) .(4,5) ,^(1,2,3) . (1,2,3,4) ,.(5) . (2,3,4,5) ,.(1) 



where L^^^^'^'* is the UV- divergent part of L'^''^^ defined by the fC-operation. Note that 

^6L?(5)' ^S?(5)' ^siS'J?)' ^£lj(i)' and are not decomposed into UV-divergent 

(2 3) 

and UV-finite parts. Note, in particular, that ' has an IR-divergent part besides a UV- 
divergent part. This is the reason why we normally avoid use of the whole L2 as a subtraction 

term and use the defined by the i^'-operation instead. For the diagram Al, however, 

(2 3) 

the IR divergences in L2 and L^llj(i) cancel each other. Thus the UV- divergence-free 
amplitude AM^i is also IR- divergence-free. The numerical integration code for the Set III(c) 
is constructed taking this observation into account. 
Substituting (1A3P in (lAip we obtain 

= AMai - LfAMsLLj, (A4) 
where Lf = L2 — L^^ is UV-finite but IR-divergent and AMg^ij = Mht,t,.t — 2Lqll{^)M2 is 



the finite g — 2 contribution from the eighth-order diagram 8LLJ 3l|, |32 |. 
Similar consideration for the diagrams Bl and CI yields 

aBi = AMbi - LfM^'^if + 2LfL,LLi3)M2 

= AMbi - L^AMslll, (A5) 

and 

aci = AMci - LfM^'^tS + 2L^Lf/2(3)^2 

= AMci - LfAMsLLK. (A6) 

From (El), (IMll . and (IMll we obtain 

c c 

J] a,i = J] AMai - LfAMsjKL, (A7) 

a=A a=A 



where AMsjkl = AMsllj + AMsllk + AMslll |31|. Note that the sum Lqll = Ull{3) + 



L%LL(A) + L^LLib) = because of the gauge invariance. 
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2. A2, B2, C2 



The diagram A2 has UV-divergent subdiagrams (1, 2, 3, 4), (2, 3, 4, 5) besides the /-/ sub- 
diagram (6, 7, 8, 9). Thus the renormahzed amphtude a^2 can be written as 

n -M /,^(l)r (2,3,4,5) . (1,2,3,4) ,.(5) 

aA2 - Ma2 - IVh ^8LLJ{2) " ^8LLG(7)^^^2 ■ l^Oj 

Diagrams 8LLJ and 8LLG are shown in Figure [31 Since A2 has no UV divergence due to 
the second-order subdiagram, we define AMa2 by 

Substituting flMj) in (IMj) we obtain 

aA2 = AM^2. (AlO) 

Similar equations hold for and 0^2- Thus we have 

c 

,2 = 5^AM,2. (All) 



c c 

a=A a=A 



3. A 5, ^5, C5 

The diagram ^43 has five forests after the /-/ subdiagrams are treated following the con- 
sideration of Sec. IIIBi Thus the renormalized amplitude 0^3 can be written as 

a A3 = Mas - ^2^^-^8LLJ(3) ~ ^i^^-^SLLH^?) ~ ^ia^'^"* -^6LL(5) 

+ M("4«'L<S,„ + Mf 'L^«'L£i>,, (A12) 

The second-order vertex renormalization constants L'^''^^ and Lg^'^^ appear in flA12p as re- 
duced diagrams, which we called implicit, and used the full renormalization constant L2 for 
them. Thus we define the finite amplitude by 

AM,, = M« - M»>LSg - mPi^^, - Af£^«>4S,„ 

+ M»>4«>iS?,5, + *ff4''*'4S,„- (A13) 

In other words, we have 

aA3 = AMa3. (A14) 

Similar relation holds for 0^3 and ac3- Thus we have 

c c 

J2aa3 = Yl ^^"3- (A15) 

a=A a=A 
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4. A4, B4, C4 



The diagram A4 has one self-energy subdiagram (2,3,4) and two vertex subdiagrams (2,3) 
and (3,4) as well as the /-/ subdiagram (6,7,8,9). Thus the renormalized amplitude aA4, is 
given by 

+ Mi'^='Br'4S,», + M<;'«M"4S,.,. (A16) 
We define the UV-finite amplitude A'Maa by 

+ M»'='Br>4a, + Mii^'^"^!^^, (A17) 

where M2* is derived from M2 by inserting a two-point vertex in the lepton line. Note that 
the renormalization constants B2 and 6m2 arising from the self-energy subdiagrams {2;a} 
and {4;d} are subtracted as a whole without breaking them up into UV-divergent and UV- 
finite parts. This is consistent only if we use the full body of the renormalization constant 
BsLLj- Otherwise, IR-singular part of B^llj and two i?2's do not cancel out each other. 
Substituting flA17p in flA16l) . we obtain 



ttAi = A'Ma4. (A18) 

The resulting 0^4 is UV-finite but IR-divergent. Separating the IR divergence of A'Ma4 
from the subdiagram {l,5;e} by the /-operation, we can write 

aA4 = AM^4 + LfiM^llf - 2LeLL(5)M2) = AM^, + L^AM^'^lf . (A19) 
Similarly we have 

asA = AMba + LfiM^'^lf - 2L,m,)M2) = AMba + L^^AM^Jf , (A20) 

and 

acA = AMba + LfiM^l'^i^ - 2L,LmM2) = AMba + L^AM^J'^J^. (A21) 
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Adding up these three results we obtain 

c c 

^ a„4 = 5Z AM„4 + L^AMsjKL, (A22) 

Q=A a=A 

noting that gauge invariance guarantees the vanishing of the sum Lqll = -^^6LL{3) + -^6LL(4) + 
LeLL{5) = 0. 

We also developed an alternative method for separating UV-divergence from Ma4, in 
which a UV-finite amplitude is defined by 

/\>'M - M ,^(1,2,5). (3,4) ,,(1,4,5). (2,3) A#(l-5) n /p\(2.3,4) ,^(1,5) . (2,3,4) 

r(1.5)D(2)/77Nr(3.4) , ,.(1,5). (2). (3,4) 



-(2,3) ^(1,5). (4). (2,3) 

"6LL(5) + ^^^2* ^6LL{5)^ 



+ M»'='Br'(B)4S',5, + A/i.'"'f 4"i£L,. (A23) 



where 



BgLLJ = BsLLj{E) + BsLLj{N), 

B2 = B2{E) + B2{N). (A24) 

The B{E) term of the wave- function renormalization constant comes from the derivative of 
the numerator of the self-energy diagram with respect to the fermion momentum p, 
while the B{N) term is the derivative of the denominator function V defined in f ill I) . For 
the second-order case, we find B2{E) = B^^ and B2{N) = B^. The relationship to the fully 
subtracted A'M^4 is thus clear and we find 

A'M^4 = A"M^4-M(''')Ai?£lf , 

^B.LLj = BsllAN) - 2B^L,LLi5) ■ (A25) 

The IR subtraction term used for A"M44 is the same one for A'M^4. As a check we evaluated 
both integrals numerically. The results are in good agreement within the uncertainty of 
VEGAS integration. 



5. A5, B5, C5 

The diagram ^45 has a self-energy subdiagram (2) and two vertex subdiagrams (1,2,3,4) 
and (4,5) besides the /-/subdiagram (6,7,8,9). The subdiagram (2) is the second-order self- 
energy diagram which contributes to the renormalization constants B2 and 61712- Taking 
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this into account we can write the renormahzed amphtude aA5 as 

+ Br>L^^,Mr + 4S,.,''4-"'. (A26) 

Applying the iC2-operation to the self-energy subdiagram (2), we obtain 

+ Br"l^,,,Mr + *,„f 4«>„iv4!-»'. (A27) 

Note that the i^2-operation yields the whole mass-renormalization constant 61712. Substitut- 
ing (1X271) in flX26D . we obtain 

= AMa5 - BfAMsLLj (A28) 
where i?^ = B2 — -B^^- Similar consideration for the diagrams B5 and C5 yields 

aB5 - ^Mb5 - ^2 ^^hlLL + ^2 ^6LL{3y^^2 + ^2 -^6LL(3)^'^2 

= AMb5 - BfAMsLLL (A29) 

_ AA/f RRA/r(l'3A5) , RRr(l,3,4) , .(5) pR.(4,5) ,.(1,3) 

ac5 - ^Mc5 - ^2 M^^^j^ + ^QLLi^y^h + ^2 ^6LL(3)^^^2 

= AMc75 - ^a^AMsLLi,. (A30) 

Adding up these results, we obtain 

c c 

^ a,5 = J] AM,5 - 52^AM8j^i. (A31) 

a=A a=A 

6. C6 

The diagram y46 has UV-divergent subdiagrams (2, 3), (1, 2, 3, 4), (2, 3, 4, 5), besides 
(6,7,8,9), and the corresponding forest structure. Thus the renormahzed amplitude a^g 
can be written as 

- M, r (2,3) ,.(1,4,5) .(1,2,3,4) ,.(5) . (2,3,4,5) ,.(1) 

aA6 - Ma6 - 1^2 M^^^j - ^g^^^(7)M2 - -Lg^^^(7)M2 

+ 4"'4S,5,Mf ' + 4^'''4ilUMf (A32) 
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Applying the i^'23-operation on Maq, we can define the UV-finite quantity AM^g as 

+ Lr'^"'4S,.,A4" + Lr'"'4'/i,M^''- (A33) 
Substituting (IA33P in (IA32p . we obtain 



= AMa6 - LfAMsLLj. (A34) 
Similar consideration for the diagrams B6 and C6 yields 

= AMb6 - LfAMsLLL, (A35) 



and 



ac6 = AMce-LfM^]:li^+LfLf^^^,^M2 + LfL^^^^^^^ 

= AMc6 - L^AMsLLK. (A36) 
From (lA34ll . ( lA35|l . and ( lA36|l . we obtain 



c 



J] a„6 = J] AM„6 - LfAMsjKL. (A37) 

7. A 7, B7, C7 

The diagram ^47 has two vertex subdiagrams (1,2,3,4) and (2,3,4,5), besides the /-/ sub- 
diagram (6,7,8,9). Thus the renormalized amplitude a^T can be written as 

_ M r (1.2,3,4) ,.(5) r (2,3,4,5) ,.(1) /.^oN 

We define the UV-finite quantity AMa7 by 

AMa-t = Ma, - ^iSfi^?^ - ^SSfi^i^ (A39) 
where Lsllg(7) is not decomposed into UV-divergent and UV-finite parts. Thus we have 

aA7 = AM AT- (A40) 

Similar relation holds for qbi and aci- Thus we have 

c c 

J]a„7 = 5ZAM«7. (A41) 

a=A a=A 
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8. A8, B8, C8 



The diagram ^48 has a self-energy subdiagram (3) and two vertex subdiagrams (1,2,3,4) 
and (2,3,4,5), besides the /-/ subdiagram (6,7,8,9). Thus, its renormahzation structure is 
similar to that of the diagram A'b: 

Applying the iCa-operation to the self-energy subdiagram (3), we obtain 

+ B^'l^i^^M? + 4?4,a4". (A43) 



Substituting flA43P in flA42l) . we obtain 

= AM^8 - E^^AM^i'Jf \ (A44) 

Applying the same consideration to the diagrams 58 and C8, and adding them to ( 1A44P . 
we obtain 

c c 

^ a,8 = J] AM«8 - ^^^AMgji^i. (A45) 

a=A a=A 

9. Sum 

Taking into account that integrals for diagrams such as Al actually represent 2x2x5 
vertex diagrams, the sum of all diagrams of Set III(c) can be written as 

C 8 

a!'"^ [Set III(c)('^'^^] = Y^Y^ Vi^M^al''^ - 3ALB2AM^y^l, (A46) 

a=A i=l 

where li refers to the open lepton line and I2 refers to the closed lepton line. ALB2 = 
Lf + Bf, and 77^ = 2 for i = 4, 7, 8, ?7i = 4 for 2 = 1,2, 3, 5, 6. 
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